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SUMMiARY 



A theoretical and experimental investigation has been 
made of the behavior of a cantilever beam in transverse 
motion with a m^ass at its tip when the root is suddenly 
brought to rest. Equations are given for determining the 
stresses, the deflections, and the accelerations that arise 
in the beam as a result of the impact. The theoretical 
equations, which have been confirmed experimient ally , reveal 
that for a beam with a given cross section and velocity at 
impact and for a given ratio of tip mass to beam mass, the 
bending stresses for a particular mode at a given percentage 
of the distance from root to tip are independent of the 
length of the beam; whereas, the shear stresses vary 
inversely with the length. 

The addition of a m.ass to the tip of a cantilever 
beam increases appreciably the stresses produced by the 
first mode of vibration but changes only slightly the 
stresses contributed by the higher modes. The tip mass 
increases the maximum bending stress much less than might 
be expected on the basis of experience with the static 
acjtion of structures. For practical engineering analysis 
tl^e maximum bending stress developed in a suddenly arrested 
cantilever beam can be found by a sim.ple addition of 
stress amplitudes in the first few modes v/ithout regard 
to phase relations between m.odes. 
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TNTROD'JCTTON 



iffhen an airplane lands, the vertical component of 
the velocity is rapidly I'odaced to rero. The shock of 
the sudden change in motion gives rise tc vibratory 
stresses in the airplane. ks a beginning in the study 
of these stresses a previous report (reference 1) dis- 
cussed in detail the behavior of a cantilever beam in 
translational motion when its root is suddenly brought 
to rest. In that paper equations are given for deter- 
mining the stresses, the deflections, and the accelera- 
tions that arise throughout the beam as a result of the 
impact. The present report extends the basic problem of 
reference 1 to include the effect of a concentrated mass 
at the tip of the cantilever beam. 

As in reference 1, the present paper is based on the 
usual engineering beam theory. In this theory the 
deflections are considered to be the result of bending 
alone, shear deflections neglected. The theory as applied 
to ordinal^ beam.s gives reasonably good results so long 
as the distance between inflection points is greater than 
a few times the depth of the beam. when this theory for 
beam action is used in vibration problem.s, such as that 
in the present paper, the results are satisfactory for 
those modes of vibration for v/hich the nodes are not too 
close together. 

This report summarizes the results of a theoretical 
solution given in appendix A and presents an experimental 
verification of these results. A numerical example for 
the calculation of the maximum, stresses near the root of 
the cantilever beam is given in appendix B. 



SYMBOLS 



E 



modulus of elasticity 



Y 



weight density of material 



coefficient' of equivalent viscous damiping of 
material 



c 



velocity of sound 
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g acceleration of gravity 

L length of beam 

I moment of inertia of cross section of beam about 

neutral axis 

A cross-sectional area of beam 

p radius of gyration of cross section of beam (^[^^ 

X coordinate along beam measured from root 

y distance from neutral axis of beam, to any fiber 

t tim.e, zero at impact 

p ooerator ~- 

o t 

n integers 1, 2, 5, etc designating a particular 

mode of vibration 

nth positive root of equation 1 + cos 9 cosh 9 
+ r9 (sinh 9 cos 9 - cosh 9 sin 9) = 0 

u 

T ratio of tip mass to beam mass — 

^ m 

Ca)„ undam.ped natural angular frequency of nth mode, 

/ 9 2\ 
radians oer second [pc 

V ^ / 

damped natural angular frequency of nth mode, 
radians per second (oo^ 



where the "frequency" is defined 




by t^n' " "n 
velocity of beam prior to impact 
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w(x,t) deflection of beam at station x and time t 



w^Jx,t) deflection of beam at station x and time t 
for the nth mode of vibration 



a(x,t) acceleration of beam, at station x and tim.e t 

a^-^(x,t) acceleration of beam at station x and time t 
for nth FiOde of vibration 

a(x,y,t) bending stress in beam at station x, distance 
from neutral axis y, and tim^e t 

av^(x,y,.t) bending stress in beam, at station x, distance 

from neutral axis j, and timio t for nth mode 
of vibration 



T(x,t) average shear stress over cross section of beam 
at station x and time t 



T^iXft) average shear stress over cross section of beam 
at station x and tim^e t for nth mode of 
vibration 

Aj^ bending-stress coefficient for nth mode of 

vibration 

shear-stress coefficient for nth mode of vibration 

C deflection coefficient for nth mode of vibration 

n 



RESULTS AND CONCLUSIONS 



Theoretical 



'J^ihen a cantilever beam with a mass at its tip is under 
uniform translation in a direction perpendicular to its 
length there is excited a theoretically infinite number 
of modes of vibration when its root is instantaneously 
brought to rest. V/ith .each successive mode, damping has 
an increasing influence upon the frequencies and amplitudes 
of vibration and, for sufficiently high modes, even changes 
the type of motion from oscillatory to nonosci llatory 
motion. In the lower modes, however, damping has little 
effect and only terms of the i^irst order in damping need 
be included in the equations. Only the equations applicable 
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to the lower modes, which alone are of importance in any 
practical case^ are presented in this section of the 
paper. For a more complete treatment of damping, see 
appendix A, 

The angular frequencies (2tt times the frequencies 
in cps) are given by the equation 



^n = 



n 



( 



where 9^ is the nth positive root of the equation 



1 + cos 0 cosh G + r9 ( sinh 9 cos G - co^h 9 sin 9)=0 ( 



In this equation r is the ratio of the tip mass to the 
mass of the beam. The values of 9^ for the first three 
modes are given in the following table for several values 
of Tl 



r 


^1 


'2 


^3 


0 


1.8751 


U.69I1I 


7 -85^8 


1 

k 


1.5738 


U.2250 


7.2815 


1 
2 


1.1+200 


I... 1105 


7-1901; 


3 

L 


1.5202 


ij..0602 


7.1539 


1 

2 


i.2i+7? 
1 . 077!!. 

. 8528 


i;.G511 
3.9826 

3.9^57 
5.9U6O 


7.1539 
7. 1026 
7.0859 
7.0802 



Figure 1 shows graphically the variation of 9^^ with the 
mass ratio r for n = 1, 2, and 5. For each value of 
n the value of 9^^, and consequently the frequency, 

decreases with increasing values of the mass ratio r. 
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Expressions for the bending stresses, shear stresses, 
deflections;: end acceJ.erations are the saiTiO as the expres- 
sions given in reference 1 for these quantities except 
that the coefficients P^, and C^-^, which characterize 

each mode, are fvinctions of an additional variable r, 
the mass ratio. The bending stress, average shear stress, 
and deflection are, respectively, for the nth mode of 
vibration; 



n 



0 p 



ii e 



2E 



S3 n (jO 



n 



(3) 



\(jO 



2 



n 



V P 
En r T- E e 



2E 



sin w 



n 



n 



2E 



sin t 



(5) 



The acceleration for the nth mode, when damping is suffi- 
ciently small, is 

a^(x,t) = w^(x,t) (6) 



The variation of the dimensi onless coefficients Aj^, 
B^, and with position along the beam. x/L is given 
in figures 2, 5, and k, respectively, for the first three 
modes, n = 1, 2, and 5 and for values of r from. 0 to 6. 
Figures 2 and 5 indicate that for all values of the mass 
ratio r the highest values of and and hence the 

highest stresses occur at the root of the beam. These 
highest or root values of A^ and Bj^ are shown for 
r = 0 and r ^ 6 in figure 5 for the first 5 modes. Root 
values of A^^ and for mass ratios between 9 and 6 are 

given in figure 6 for the first 5 modes. Both figures 5 
and 6 show that the addition of a mass at the tip of the 
beam (r 0) increases appreciably the values of the stress 
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coefficients and Bj.;^ for the first mode (n = 1) but 

has a very small effect upon these coefficients i'or the 
second and higher m.odes. 

The tip mass increases the maxim\im bending stress 
much less than might be expected on the basis of static 
considerations. For exam^ple the addition of a tip mass 
6 times the m.ass of the beam, increases the mass moment 
about the root 1200 percent whereas the first mode 
bending stress coefficient A-]_ is increased only l8J,' per- 
cent (from 1.566 to U.LSO). (See fig. 5.) 

The maximum values with respect to time of o^(x,y,t) 
and T"^(x,t) associated with the nth mode of vibration, 
when the effects of damping are neglected, are 



Equations (5) and (I4.) for bending and shear stress, 
from which equations (7) and (6) are obtained, and equa- 
tions (5) and (6) for deflections and accelerations give 
the values associated with the nth mode of vibration. 
Since all modes of vibration occur simultaneously the net 
results are the superposition of the effects of all modes. 
This superposition gives the following equations: 

/ X.(jo 



a(x,y,t) = - E ^Aj e ^ sin co-j^ t 



a^(x,y) - A^ ^ ^ E 



(7) 



V p 
- - E 

c L 



(8) 



2 




t 




e 



2E 



t + 



(9) 



2 



8 
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For average shear stress, 

/ A,00-, 

~(x,t) = X i 3 e sin t 



Xa>2 \ 



+ e 2^ sin CO., t . . . / (10) 



For deflection. 



w(x,t) = ^ ^ ® ^1 ^ 

t 



+ e sin 00^ t + • . .y (11) 



For acceleration, when damping is sufficiently small, 

2 



V l2 I. ^ 



a(x,t) =--^ — ^i^i^e 2E 3-^^ 00^ t 



4- C2 ^2 ^ (JO2 t -f . . . / (12) 



For a beam v/ith a given cross section and velocity 
at Impact, the equations for bending stress reveal that 
at a given percentage of the distance from, root to tip 
and for a given mass ratio, the bending stress for a 
particular mode is independent of the length of the beam. 
The equations for shear stress reveal that the shear 
stress at any station varies inversely with the length 
of the beam.. 
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Experi mental 

The apparatus which was used to provide for the 
instantaneous arrest, of a cantilever beam is shown in 
figure 7* th5 s apparatus two cantilever beams are 

formed by centrally clam.ping a steel tube In a heavy split 
block. The block is attached to a carriage which is 
permitted to run with known velocity over a horizontal 
track. The carriage is accelerated by a weight acted 
upon by gravity in the initial portion of the run and is 
kept in a state of uniform translation by an additional 
small weight used to overcome friction in the latter 
portion of the run. Instantaneous arrest is achieved 
by permitting a tapered plug projecting from the 
carriage to ram into a fixed chuck. The effect of a tip 
mass v;as studied by Increasing the weights on the tips 
of the beams in successive tests. The velocity at Impact 
and the dim.ensions of the cantilever beam:S are given in 
appendix B. 

The apparatus described herein provides for a miuch 
more rigid clamping of the tube and gives a better control 
over the instantaneous arrest than the apparatus described 
in reference 1. With the more rigid clamping, less 
oscillatory energy was lost by the cantilever beams to 
adjacent parts of the apparatus. The damping present, 
therefore, m.ore nearly approached the dam.ping of the 
miaterial of the beam. 

Extreme fiber bending stresses near the root of each 
cantilever beam were measured by means of electrical 
strain gages and a recording oscillograph as described 
in reference 1. A typical record of the bending strains 
at the roots of the two cantilever beams with the mass 

ratio r = p is shown in figure 8. No measurements were 
k 

made of the shear stresses since their values were too 
small to be m.easured accurately in the presence of the 
vibrations set up by the rolling of the carriage. 

The three quantities that were obtained from the 
tests - the frequencies of the first three modes, the 
maximum contribution of the first mode to the total 
extreme fiber bending stress at the root, and the 
maximum extreme fiber bending stress at the root - 
are plotted against the mass ratio r in figures 9 
and 10 for comiparison with the theoretically computed 
values. Since inherent local variations in the beam 
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oroperties do not appreciably affect the frequencies, 
which are associated with the over-all action of the 
beain, and since frequencies are easy to measure, it is 
reasonable to expect the observed good agreement between 
theoretical and experimental frequencies, (See fig. 9.) 
^Ah.en consideration is given to the fact that stresses are 
directly affected by the local variations in the beam 
properties and are not readily susceptible to instantaneous 
accurate measurement the observed agreem.ent between the 
experimental and theoretical stresses is also considered 
to be satisfactory. (See fig. 10.) 

The coiit ribution of the first mode to the total 
stress vvas estimated from the records. (See fig. G. ) 
It is clear from figare 10 that the first mode contributes 
more than half of the total stress. It is also clear from 
figure 10 that for practical engineering analysis the 
maximum bending stress developed in a suddenly arrested 
cantilever beam can be found by a simple addition of 
stress amplitudes in the first few miodes (in this case 5) 
without regard to phase relations between the modes. 



Langley Kem.orial Aeronautical Laboratory 

National Advisory oommlttee for Aeronautics 
Langley Field, Va. , November $0, 19l\h 
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APPE^IDIX A 



TFEORETIC'^L DERIVATION 



Ge n e r a 1 an a 1 y s Is . - Consider In equilibrium a uniform 
cantilever beam with a mass at its free end. If the root 
of the beam is sr.adenlj/ disturbed, as by a shock, in a 
direction perpendicular to its length, the beam will be 
set into damped bending osci llat:^ ons • The equation of 
motion for these bending oscillations is given by the 
differential equation (reference 1) 



With the use of the notation = ^~ equation ( Al ) may 

be written 

This partial differential equation is reduced to an 
ordinary differential equation of the kth order by 

writing p — - , thus. 



1 + ~) ~T + w - 0 (A3) 



The general solution of equation ( A5 ) is 

w = p cosh G-J + Q sinh e~ + R cos 6- + S sin 67 (aL) 
L T, L ij 
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v;here 




The coefficients P, Q, R, and S are to iDe determined 
from the boundar;^/ conditions. The case under considera- 
tion is that of a cantilever with a mass at its free end 
moving ^ivith uniform velocity- v and having its root 
brought instantaneour3ly to rest. The boundary conditions 
for this case are 




The fourth boundary condition, which is an applica- 
tion of Newton^ s third law, equates the shear force at 
the tip of the beam to the inertia force of the tip mass. 
The velocity of the root as given by the first boundary 
condition is represented graphically in figure 11. 
Following the procedure adopted in reference 1, the 
solution \r^ill be obtained for the boundary condition 



and to the resulting velocity will be added the constant 
velocity v. 
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With the aoplication of the boundary conditions to 
equation (aIi), the operational solution for the velocity 
(that Induced by -vi ) Is found to be 

pw = 

^ 1 + cos 9 cosh e + r0 ( sinh e cos 6 -cosh 6 sine) 

(A5 

where 

^(^r) " ^ ^^^^ e)(^osh 0^ + cos e^^ 

+ sin •9' sinh G Tcosh 0 — - cos 0--^ 

+ ( sinh 9 cos 0 + cosh Q sin 0) ^sin 0~ - sinh 0~ ^ 



+ 2r0 



X -yr 

sinh 9 cos 0 cosh 0^ - cosh 0 sin 0 cos 9^ 



/ X X N 

+ cos 0 cosh 0(^sin 0~ - sinh J 



and r is the ratio ~= _ Tip ma ss ^ interpretation of 

Beam mass 

this operational expression and addition of the con- 
stant velocity v gives for the total velocity 



- ^ t / n 
6w(x,t) o V" ^^ 2E / \ 



(a6)'^- 
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where 

Qj.^ nth positive root of the frequency equation 

1 + cos 9 3ndiQ -f r9 ( slnh 9 cos 6 - cosh 0 sin 0 ) - 0 
(all roots, namely ±9^, *^^n^ have been considered 
in the interpretation) 



^n " ~ 



undamped natural frequency of nth mode, 
radians/sec 



2 2 



X CO 



n 



damper^ natural angular frequency 
2 of nth m.ode, radians/sec 



cos 0^ cosh 0jT^)(cosh 6^ 

r) 



I, ^ 9^ --j 



- COS 9n 



! + sin 9^ sinh 9^ t cosh ^ 
^ ' 3inh 9^^ cos 9^ 

^ cosh 9^ sin 9n 



( sinh 9^-^ 
^ 2r9n 



9^ sinh 9^ 

+ cosh 9n sin 9^) (sin ^ 



.J. sinh 9n ~) 



sinh cos 9^ 



/ . + cos Or, CO 

^ e^[(L-r)(sinh8nCos 6^ 



cosh 9n Y" 
cosh 9j^ (sin 9jf - sinh 9^, j- j 

>-2r9^ sinh 9n 



- cosh sin 9;^)- 



> 

in cos Sn^ 
sin gjj 



Integration of equation (a6) vlth respect to time vath 
the conditloTi (w)^_q = 0 gives for the deflection 

when t > 0 



, t) = 2v\ -^E sin03„'ti 

^n' 



w( X 

n^ 



n=l 
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where 

2P 

C 



n 2 

The contrxbutlon of the nth mode to the deflection 

is 

2 ^^n^ 

V L ^ — t 

^n(^'t) = - -- -p._J--._- e 2E si^c,^. ti (a8) 

/I P_ 

Vi/hen — > 1 equation (a8) may be put in the form 

2 ^^t 
w^(x,t) =^^C^ L e 2E 3i^^,, ,ti (A9) 

~ - 1 



where now 



= ^nl/ ~ - 1 



The form indicated by equation (a8), where ~^ < 1 is 

2ji ' 

characteristic of the lower modes and represents dam^ped 
oscillatory motion. The form indicated by equation (A9), 
where --ii > 1 (damping greater than critical), is 

characteristic of the higher modes and represents sub- 
sidence m.otion. 

FroiT! equation (a6) for the velocity and equation (A?) 
lor the deflection, the complete behavior of the cantilever 
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may be determined. The quantities of interevSt are tbe 
bending stresses, the shear stresses, and to some extent 
the accelerations. ^ATiere damping is present;, the equa^ 
"cions representing the contribution of the nth mode to 
these quantities may be given in the tv/o forms indicated 
by equations (a8) and ( A9 ) • subsequent equations, 

hov^ever^ only the form indicated by (Ao) is given, because 
it is characteristic of the miodes v/hich are of practical 
importance. 

Bonding stresse s . - The bending stresses a(x,y,t) 
at any fiber a distance y from the neutral axis are 

a(x,y,t) = Ey = E ~ ^ ^n'*^ 

6x'"" n--l 

where 

'"(1 + cos 9^ ^^osh 9 ^) ^cosh 9n^ " ^"^^^n-r} 
+ sin 9^ sinh 9^^osh 9^ ^ cos 9^^ J 



n 



> -( sinh 9^^ COG + cosh 9 sin 9^) Tsinh 9^ ^ ^n ^) 
+ 2r9n jsinh 9n cos 9^ ^osh 9^-—- + cosh sin 9^^ cos 9^^ ^ 
- cos 0^ cosh 9n(sinh 0^^ f + sin 9n y 



^ r 1 

8^1 (1 + r)( sinh 9^^ cos 9^ - ''^osh 9^ sin 9^) ~2r9n sin sinhgj 

The contribution to bending stress of the nth mode is 
an(x,y,t) = E J ^ An e co^ ' t i 



Average shear stresses . - The average shear stress 
on the cross section is 
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n ^--^ — t 



V 



sin ^3;^ ^ t i 



where 



< 



(1 + cos cosh 9,,)(siBh 8^ 2: 0^ £j 

sin sinh 9^(sinh 0^ f - sin 0^ 



-(sinh cos 9„ + cosh 9^, sin 9n)(cosh 6^ - + cos 8^ ^) 
+ 2r9^rsinh cos 8,^ sinh 8^ ^ - cosh 9^ sin 6^ sin Q^^ 
- cos 8^ cosh 9 



> 



v.- uu. o -usn y^(cosh f + cos gj] ^ 
(1 + r) (sinh 9„ cos 9^ - cosh 0^^ sin Q^) - 2rQ^ sin 8n sinh 8^^ 



mode 3^s^ ^contribution to average shear stress of the nth 



\(jO 



1 £ 
c L 



n 



^ - - e 2E 



;in 00^^ ' t J 



Acceleraci_ons.- From equation (a6), v/ith the aid of 
the relation 

pP(t)i = F(0)pi + F'(t)i 
the acceleration anywhere on the beam Is found to be 
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■6t2 



c P 



r=l ^ J 



n 



.op ^ 



! s 1 n '-^3^ ^ t 



n 



11 - 



_i i 

2E'~ 



COP, ' t j i 



/ 



With the aid of "che orthogonal properties of the func- 
tlons i; j is possible to show that the quantity 



00 

2 ^ -^(^n ^) ~ ^^^'^^^^ to zero 



when 0 < T 1 1. 



n=l 



At ^ equal to zero the quantity 2 > kn ^) Is zero 

n=l 

and only the term -vpi remains. This term indicates 
that at t = 0 there is at the root an infinite accelera 
ti on of zero duration. 

The contribution to acceleration of the nth mode is 



r2 1 - — ^- 

c p "^n r \2<j) 2 ® 



'1 - 



a 



sin ' t 
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Comparison with the exDression for w,^(x,t) (equation (a8)) 
shows that the acceleration of each mode is out of phase 
with the deflection. V.hen dam.ping is sufficiently small, 
however, the relation between the acceleration and the 
deflection reduces to the well knov/n result for undamped 
vibration 
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APPENDIX B 



N^imeiloal Example 

Problem. - To calculate the theoretical bending and 
shear stress at the position of strain gages on the steel 
tube beam used in the experimental investigation of the 
present papei* for a mass at the tirj equal to the mass of 
the beam, r = 1, 



Length of beam., L, In 29.8I4. 

Outside diameter of tube, in„ • . . • , l.OOJ 

Distance to extreme fiber, maximum value 

of y, In 0.502 

Wall thickness of tube, in 0.029 

Radius of gyration of cross section, p, in. . . . 0.5^+5 
Distance fromi root of beam to strajn 

gages , X, in '^•50; 

Modulus of elasticity, E ( ass-omed), psi .... 29 x 10; 

Velocity at im.pact , v, fps ^1-7^ 

Velocity of sound in steel, c, fp^. l6,600 



The effects of damping will be neglected so that equa- 
tions (7) and (8) may be used. From the foregoing data 
of this problem 



11.^ = jj.500 DSi 
CP 



- 7 S = 5p-b psi 
c L 

The values of A.^ and Bj^ for the different modes are 

obtained from, figures 2 and 5 for r = 1 and 7 - 0.0167. 

The computed stresses are given in the follo7;ing table 
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Mode , n 


An 


— — 

Bn 


a 

(psx ) 


(psi ) 








(ii^OO X An) 




1 


^.13 


2.I4.O 


9800 


86 


2 


.89 




liOlO 




5 


.U3 




2160 


li;5 


Sum. of f'irst 


three 


15,970 


373 



stress amplitudes 



An approximation of the maximu-in total stresses can be 
obtained by adding the stress amplitudes for the first 
several modes as indicated. 
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Figure I." Variation of 0^ wi+h mass ratio r 

for n " 1,2, and 3 • 
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Figure 2.« Variation of bending stress 

coefficient Ap with )/L. '^"^ 
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Figure 3. ~ Variation of shear stress 
coeffic»en+ Bp with x/L . 



NACA MR No. L4K30 




NACA MR No. L4K30 



at 

root 



2 - 



2 3 

Mode , n 



JEL 



NATIONAL ADVISORY 
COMMIHEE FOR AERONAUTICS 



at 
root 



3 
2 



2. 3 

Mode , n 



r = o 



0 



r= 6 



Figure 5. -Values of bending - stress coefficient A^, 
ond shear stress coefficient Bn Q"*" ''oot of 
cantilever beam for moss ratios r * O and 6. 
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Figure 6.- Varia+ion of bending - stress 
coefficients Ap and shear- stress coefficients Bn 
at root of cantilever beom with mass rotio r 
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Figure 9 Comparison of theoretical with 

experimental frequencies as mass ratio r increases. 



3b 



NACA MR NO. L4K30 



16 




Modes 1 + 2 + 3 

Addition of stress 
I amplitude 

+ 



Addition of stresses with 
proper regard to phose 



4 




X Maximum recorded stress 

G riQximum stress in mode I 

NATIONAL ^DVISORV 
COMMI TEE fOF AERONAJTICS 



.8 



1.2 



1.6 



2.0 



Figure 10.- Comparison of theoretical with 
experimental root bending stresses os 
mass ro+io r mcreases . 
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Figure II.- Graphical representation of the 
discontinuous function v-v/. 
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